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We carry out a calculation of the phonon-drag contribution Sg to the thermoelectric power 
of bulk semiconductors and quantum well structures for the first time using the balance equation 
transport theory extended to the weakly nonuniform systems. Introducing wavevector and phonon- 
mode dependent relaxation times due to phonon-phonon interactions, the formula obtained can be 
used not only at low temperatures where the phonon mean free path is determined by boundary 
scattering, but also at high temperatures. In the linear transport limit, Sg is equivalent to the 
result obtained from the Boltzmann equation with a relaxation time approximation. The theory is 
applied to experiments and agreement is found between the theoretical predictions and experimental 
results. The role of hot-electron effects in Sg is discussed. The importance of the contribution of Sg 
to thermoelectric power in the hot-electron transport condition is emphasized. 
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'nI" ■ I. INTRODUCTION 

^ ■ 

, Recently the study of thermoelectric power has attracted intengified interest both theoreticaUy and experimentaUy. 
T-H in various systems: mesosflopic quantum dotad'EI, quantum wiresu, heterojunctions and quantum well structuresQ Ell 
as well as bulk materialstirEJ. Almost all of the earlier theoretical treatments were based-pp-.ijlie Boltzmann equation 
0^ in analyzing the-tjUfEusion component ofj-thermoelectric power in macroscopic systemstjEO'cl and the phonon-drag 
component of itll2l3ia. In a recent lettertil Lei pointed out that thermoelectric power can be conveniently analysed 
entirely within the framework of balance equation approach. He discussed thermoelectric power of both bulk materials 
and quantum wells in the piEesence of carrier heating with a strong applied electric field using balance equations for 
weakly nonuniform systemalj. His calculation indicates that hot-electron effects on thermoelectric power may not 
' only change its magnitude, but_aJso may change its sign at higher electric field. His result was confipned very recently 
Q ] by Xing, Liu, Dong, and Wan^ usiug the non-equilibrium statistical operator method of ZubarevEJ jointly with the 
O ■ Lei-Ting balance equation approachE3. However in assuming that the phonon-drag contribution tp_.thermoelectric 
^ ' power may he jieglected in the range of electron temperatures of interest for hot-electron transportllj, both of their 
• ^ • treatmentstilOr-Hiissed this contribution, which is [kBOwp to be very important in linear transport in both bulk 
' semiconductorsEJ and two-dimensional (2D) systemad uM. It is usually so large that the diffupion contribution would 
not be discernible unless measurements were extended to very high or very low temperaturesc3. 

In our opinion this assumption is questionable. The analysis reported here includes the phonon-drag component 
Sg of thermoelectric power using the balance equation approach, which has the advantage of easy inclusion of hot 
electron effects. In this, we clarify the importance of the contribution of Sg to thermoelectric power in the hot- 
electron transport condition. Our consideration are applicable to the regime in which the electron drift velocity is 
lower than the sound velocities, for materials having high impurity concerntratinos-and/or intermediate electric field 
strength. In fact, our result shows that contrary to the assumption of Xing et aLl3, Sg is markedhi, enhanced at low 
lattice temperature under the indicated conditions. As the calculations of LeiM and Xing et altB both show that 
the hot-electron effect does not change the absolute value of the diffusion component of thermoelectric power very 
much, it is therefore clearly important to include Sg in the hot-electron transport condition, provided that the lattice 
temperature is not too high. 

This paper is organized as follows. In Sec. II, we present hydrodynamic balance equations with the phonon-drag 
effect included. In Sec. Ill, we compare our theoretical prediction of Sg with experiments in the linear transport limit. 
Finally, in Sec. IV, we discuss hot-electron effects on Sg and present our conclusions. 
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II. HYDRODYNAMIC BALANCE EQUATIONS INCLUDING PHONON-DRAG 



The balance equations for electron transport undetthe influence of an electric field E and in the presence of a small 
lattice temperature gradient VT, can be written asEj: 

^+V-(vn)=0, (1) 

+ vVv = -- + — + - 2 

Ot 6 mn m mn 

^ + V • Vu = -^^(^ • V) - W - V (f, + fep) . (3) 

In these equations the carrier drift velocity v, the electron temperature Tg, the average relative electron energy u, 
the carrier density n and the chemical potential /i are all weakly dependent on spatial coordinates due to the small 
lattice temperature gradient, so that their spatial gradients are all small. Eqs. (|l|)-(§) are accurate to the first order 
in these small quantities (all the second and higher orders of spatial derivatives are neglected). The frictional forces 
due to impurity scattering and due to phonon scattering fgp, and the energy transfer rate from the electron system 
to the phonon sjptem w, are all functions of local quantities n, v, and T. fi has exactly the same expression as in 
the uniform casalil. i^p and w can be written as {h — kg = 1) 

fep = -47r^ |Af(q, A)pq(5(ek+q - £k + - q • v) 

kqA 

X {nqA[/(ek,re) - /(ek+q,T,)] - /(ek,re)[l - /(ek+q,re)]} , (4) 
w ^ 47r^ |M(q, A)pf7qA (5(£k+q - £k + f^qA - q • v) 

kqA 

X Mfiek,Te) - /(ek+q,Te)] - /(£k,Te)[l - fie^+^,T,)]} , (5) 

where £k = k'^ /2m is the energy dispersion of an electron for an isotropic parabolic band system, flq\ denotes the 
energy of a phonon of wavevector q and branch A, and /(sk, 7e) = l/[exp(ek — M)/2^e + l] is the Fermi distribution func- 
tion at the electron temperature, n^x is the occupation number (distribution function) of phonons in mode qA, which 
is weakly nonequilibrium due to the existence of the small thermal gradient. It is just this gradient that leads to a net 
flow of phonons in one preferred direction, which then imparts momentum to electrons via electron-phonon scattering. 
Therefore, in the balance equation approach, the phonon-drag contribution to thermoelectric power is embedded in 
the nonequilibrium phonon distribution function nq\ and thereby in the electron-phonon-scattering frictional force 
fep, which was previously taken to be the equilibrium Bose distribution function n{ilqx/T) Lil[exp(fiq>/r) — 1] in 
Refs. |l2| and |ri|. M(q, A) stands for the dynamically screened electron-phonon matrix elementEilc3. Considering the 
relation 



'5(ek+q - £k + r!qA - q • v)[/(ek, Te) - /(ek+q, T,)] „( "qA q-v ^ ^ 



^(£k+q - ek + f^qA - q ■ v)/(ek, T,)[l - /(ek+q, T,)] 

Eqs. (^ and (||) can be rewritten as 

fe, = -2 5] q|A?(q, A)pn2(q, l^q^ - q • v)[nqx - n(^^3^— ^)] , (7) 

«; = 2^r!qA|Af(q,A)|2n2(q,r!qA-q-v)[nqA-n(^^3^—i:^)] , (8) 

qA 

in which 112 (q, ^^qA — q ' v) stands for the imaginary part of the density-density correlation function 

n(q,^) = 2yi^i^i^±2l^^l^/(filll^. (9) 
^ ' ^ £k+q - ek + w + iO+ ^ ' 

The phonon distribution function n^x for the nonuniform system can be determined as follows. Its total rate of 
change can be written as 
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+f%i-) +f%| , (10) 



dt ^ \ at J, \ at y„„ V at 

where the diffusion term is 



-Vp(q,A) • VnqA . (11) 



V dt J, 

To the lowest order in the thermal gradient, the distribution function nqx on the right hand side of Eq. ( [Til ) can be 
replaced by the equilibrium Bose distribution function: 

-v,(q,A).Vn(%^). (12) 



V dt ' ' T 

In these equations the dependence of n^x on spatial coordinate r has been suppressed and is understood; Vp(q, A) 
stands for the group velocity of a phonon in mode qA, and is given by VqfiqA- The second term on the right hand side 
of (|l^) is the rate of change of n^^x due to clcctron-phonon interaction, which dpscribes the absorption and emission 
of phonons by electrons in transport. It has been proved that this term readf£3c3 

=2|M(q,A)pn2(q,l]q.-q.v)[nq.-n(^^3^^^)]. (13) 

This effect, together with the diffusion of phonons, tends to drive the phonon system out of equilibrium. On the other 
hand, there is a trend to drive these nonequilibrium phonons back towards equilibrium, which is described by the last 
term in Eq. ([Tc|): 

V cit Tp{q,X) T 

Here Tp(q, A) is a wavevector and mode-dependent relaxation time, which is assumed to be determined by (a) boundary 
scattering 

1/tsa = v,x/L , (15) 

(where Wpi_denotes the speed of sound in mode /L-and L is of the order of the macroscopic dimension of the 
specimenotEl) , and by (b) phonon-phonon processescZltj: 

l/Tpp(q,A) = AAT3f]2^, (16) 

whence 

l/rp(q, A) = 1/tbx + l/rpp(q, A) . (17) 

In the steady state, the occupation number for each phonon mode is constant, ie. dn^x/dt = 0, and Eqs. (p^)-(p^) 
lead to 

( anq^y , 1 „r"q^'i -I 1 / "q>-q v \ 

V dt T^{ci,X)^^ T ^ + r(qA,nq;,-q.v)"t- > 
"qA = 11 : (1C5) 

■r"p(q,A) T(qA,OciA-q-v) 

in which we have introduced a phonon-mode dependent and electron-drift- velocity dependent inverse electron-phonon 
scattering time 

1 -2|M(q,A)pn2(q,l}qA-q-v) , (19) 



T(qA, fiqA - q • V 



and the phonon drift component ^ ^^^^ ^ has been replaced by Eq. (p^). Substituting the distribution function (|T^ 
into Eqs. (Q) and and taking the definition ( [l9| ) into account, we have 



fep = fg + fp , w = Wp , (20) 
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where fg is the frictional force due to phonon drag; tp and Wp are the frictional force and energy transfer rate due to 
the phonon absorption and emission. They are defined by 



qA 



= E q^^Y^TTo rt'^^^) - "(^^^^^)] ' (22) 

^ rp(q, A) +T(qA, SiqA -q-v) T Te 

qA 

»-'"'-5:«,^^STjT7sW^l"'%i)-»(23^^ (23, 

qA 

We consider a small lattice temperature gradient along the x direction: VT — (Va;T, 0,0). Then, in addition to 
the drift velocity (current flow) and the applied external electric field along the y direction, there will be a small drift 
velocity (current) and a small electric field in the x direction: v = {vx,Vy,0) and E — (E^, Ey,Q). Here, and 
are small and [the spatial variations of all quantities are assumed to be along the x direction only. Following the same 
scheme as Leilj by considering the particle number, force and energy balance equations (|l|)-(g|) to first order in the 
small quantities, we obtain, in the steady state 

--^V^u+^ + ^ + ^^O, (24) 
6mn mn m ran 

neEy + fy=0, (25) 
w + vyfy=0, (26) 

which state the force balances in the x and y directions, together with the energy balance, respectively, and f = f ; +fp. 
For small Vx, fx is proportional to it, and p = —fx/in^s'^Vx) is the resistivity in the x direction in the presence of 
drift velocity Vy in the y direction. Then Eq. ( p^ may be rewritten as 

J.^^-l^ + ^. (27) 
p 6 nep nep 

Noting that the energy u and number n densities are given separately by 

u^2j2e^fi'-^) , n^2E/(^), (28) 

k ^ k ^ 

we can rewrite Eq. (^^ in the form 

jx = L^\Ex ~ Vxfi/e) + Lf{~VxTe) + Lf{~VxT) , (29) 

in which L^^ = 1/n. L]? = — V'J^'^'^/'l — CI denotes the coefficient of the electron diffusion contribution to the current 

with C = p/Te and Fy{y) — a;''[exp(x — y) + l]^^dx. L]^ is the coefficient of the phonon drag contribution which 
is obtained from Eq. (21) as 

= ^^^^OjA) (30) 

9 nep^g2 ^2 ^p-i(q,A)+r-i(qA,l^qA-c.o) t'' 

where we have used the Debye-type spectrum for acoustic phonons: fJqA = Vgxq with Vs\ being the longitudinal sound 
speed in A mode and lhq f^fnVy The thermoelectric power S of hot electrons in the presence of a strong current flow 
in the y direction is thenll3'0 

ri2 xrp ^12 

The diffusion component Sd has already been obtained by Leiii and Xing et a/.lli. When there is no current flow in 
the y direction {ST^/ST = 1), it matches the result obtained from the Boltzmann equation with a constant relaxation 
time assumption in Ref. 0. Furthermore, 
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is the phonon-drag component of thermoelectric power. Thus, Eqs. (|3]1)-(^2|) yield the thermoelectric power directly. 

The hydrodynamic balance equations (|l])-(^) for weakly nonuniform systems are also applicable to 2D systems 
(quasi-2D elections coupled to bulk phonons) if the prefactor | in Eq. (H) and | in Eq. (^) are replaced by unity and 
2, respectively^. Moreover |M(q, A)|^ in Eqs. (Q) and (||) is replaced byiM(q, A)p|/(«g^)p. Here, I{iqz) stands for a 
form factor associated with confinement of the sample in the z directionEiJ. Specifically, for a quantum-well structure 
it reads 

\I{tq.)\'^7T^sm^y/[y\y^-7rY] (33) 

with y = qza/2 and a is the width of the quantum well. Thus, the diffusion component of thermoelectric power Sd 
for a 2D system takes the form 



e 



2^^i(C) 



FoiO 



The phonon-drag component Sg shares the same expression as that of Eq. (|32|). We may put this equation in a more 
symmetrical form when there is no current flow in y direction, and correspondingly luq = 0: 

S 1 y gji ^qA 1 r/f^'^^l (35) 

qA ^ Tp(q,A) 

with q = {(\\\,qz) = iqx,qy,qz)- X^is equation can be shown to be equivalent to that derived from the Boltzmann 
equation by Cantrell and ButchenJ if "1" in the denominator of Eq. ( ^ ) is omitted. This correspondence arises 
because Cantrell and Butcher use an approximation (ie. see Eq. (41) in the first paper of Ref. ^ which is analogous 
to neglecting the electron-phonon scattering time l/r(qA, fiq^ — q • v) in our formula (p^. If we eliminate this 
approximation in their derivation, then their final result may be proved to be equivalent to Eq. (|35|). In our opinion, 
this approximation is valid only at low temperature. If temperature is high, this approximation yields a result that is 
too large. 

III. COMPARISON WITH EXPERIMENTS IN THE LINEAR TRANSPORT LIMIT 

In this section we apply the present theory to GaAs/GaAlAs quantum wells. Both longitudinal and transverse 
acoustic phonon modes are coupled to electrons via deformation potential and piezoelectric interactions. Only the 
longitudinal mode gives rise to deformation potential coupling, whereas both longitudinal and (two branches of) trans- 
verse modes contribute to piezoelectric interaction. Their matrix elements are given in Ref. The GaAs/GaAlAs 
material parameters used in our calculations are as follows: volume density d = 5.31 g/cm'^, effective mass m = O.OTme, 
transverse sound velocity Vgt = 2.48 x 10^ m/s, longitudinal sound velocity Vsi = 5.29 x 10'^ m/s, acoustic deformation 
potential S = 8.5 eV, piezoelectric constant ei4 = 1.41 x 10^ V/m and low-frequency dielectric constant k = 12.9. We 
note that some of the values of these parameters are still uncertain ( eg. in the literature the values for S range from 
7 to 16 eV). Nevertheless the above listed values have been used frequently in balance jeauation theory in a variety of 
transport problems and the results obtained are in good agreement with experiments. The coefficients A\ in Eq. 
( |l6| ) can be obtained by fitting to experimental data, and a particular choice will be identified below. j-. 

In Fig. 1 we first fit Eqs. (31), (p3) and (^Sf) to the experimental data of Sample 1 of Fletcher et ala, where 



the electron sheet density is = 1.78 x 10 m~ , and the length between the contacts used to determine the 
thermoelectric power is L = 0.98 mm. As there is no external electric field in the experiment, 5Te/6T = 1, ujq — 
and Te = T. In our fitting we take Ai = 1.68 x 10~^^ s K~'^ for longitudinal phonons and At = 4.36 x 10~^^ s 
for transverse modes in Eq. (p^). For comparison purposes, the thermoelectric power components due to deformation 
potential coupling Sgi , those due to longitudinal piezoelectric interaction Sgpi and those due to transverse piezoelectric 
interaction Sgpt are also plotted in the same figure. We then apply the same parameters to the calculation of S for 
the other two samples of Fletcher et alB. The plot of S vs. T in Fig. 2 is for electron sheet density N — 4.8 x 10^^ 
m~^, with length between the contacts L = 0.63 mm, and in Fig. 3 it is for N = 6.82 x 10^^ m~^, with L = 2.98 
mm. The widths of the quantum wells in Fig. 1-3 are aU 100 A. From the fi gures we can see that our theory is 
in good qualitative agreement with the experimental data. Furthermore, our results in Fig. 1-3 show that, at low 
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temperature, the two branches of transverse phonons dominate the drag effect through piezoelectric interaction with 
electrons, whereas the longitudinal phonons contribute to Sg mainly through the deformation potential coupling. It 
is also seen that the phonon-drag component dominateSjihermoelectric power in the experiments. 

We also apply our theory to a very recent experimentEJ, where the electron density is = 3.6 x 10^^ m^^, which 
is of an order of magnitude larger than that of the previous samples. The width of the quantum well is a = 100 Aand 
L = 2.5 mm, as determined in Ref. |l^ using experimental data for thermal conductivity. We use the same material 
parameters and coejjififnts Ax as in the previous cases. Figure 4 exhibits a comparison of the experimental data with 
our theoretical resultcEI. Again the theory fits well qualitatively with experiment. However, this time the theoretical 
prediction is smaller than the experimental data even at low temperature. _ 

This experiment has also been explained by Fletcher et al. in the same paperE2l using a formula derived from the 
Boltzmann equation by Cantrell and Butcheiu. The material parameters they used are chosen to fit the experimental 
data at low temperature. However these parameters faiLto fit the earlier experimental results shown in Figs. 1-3 
at low temperature. In the calculation of Fletcher et al£3, the relaxation time of phonons is treated as a constant 
and is determined from the thermal conductivity measured in their experiment. Even though the temperature in 
the experiment ranges from low temperature to very high values, they continue to adopted the approximation of 
Cantrell and ButcheiCI throughout, whereas we believe it to be valid only at low temperature. In Fig. 5 we compare 
our theory with that of Fletcher et al., wherein the chain curve is calculated from our theory and the dotted curve is 
the theoretical prediction of Fletcher et al.. It is seen that the peaks predicted by the two theories occur at almost 
the same temperature. Nevertheless, the approximation used by Fletcher et al. yields the height of the peak of their 
theory to be larger than it should be, because this approximation fails at high temperatures. This may be understood 
from the fact that if we include the same approximation in our formula {ie. omitting "1" in the denominator of Eq. 
(p5|)), the result obtained {ie. solid curve) predicts the same peak height found by Fletcher et al.. 



IV. HOT-ELECTRON EFFECTS ON THE PHONON-DRAG COMPONENT Sc AND CONCLUSIONS 



This theory may be easily used to calculate thermoelectric power in the presence of a strong external electric field. 
In fact, althpugh the diffusion component Sd may be negative within a small low-lattice-temperature range at higher 
electric fieldtil, the phonon-drag component Sg is still positive. Moreover, it should be even more important to include 
the phonon-drag contribution to thermoelectric power in the hot-electron transport regime. Actually, Sg is markedly 
enhanced at low lattice temperature as shown in Fig. 6, for the sample N — 1.78 x 10^^ m^^ and in Fig. 7 for the 
much higher electron density sample N = 3.6 x 10^^ m^^. Our calculated results shown in Fig. 6 and 7 assume high 
impurity density, such that the electron drift velocity in hot electron transport is smaller than the sound velocities 
{ie. we discard in Eq. (32)). In Fig. 6, Sg is plotted against lattice temperature T as a solid curve for hot-electron 
temperature T,. = T -\- 200 K, as a dashed curve for T,. = T + 100 K and as a chain curve for Tg = T -I- 10 K, with 
all the other material parameters being the same as those used in Fig. 1. For comparison, Sg in the linear transport 
limit is also plotted in the same figure (Fig. 6) as a dotted curve, with part of it obtained in Fig. 1. From the figure 
one can see that Sg is markedly enhanced by hot-electron effects in the low lattice-temperature regime (LLTR) with 
T < 20 K. However, in the high lattice-temperature regime (HLTR) with T > 20 K, is slowly reduced as Tg 
increases. Another feature to be noted is that for T < 6 K and T > 10 K, Sg, for = T + 200 K (solid curve), is 
a bit smaller than it is for = T + 100 K (dashed one). These properties can be well understood from the inset of 
Fig. 6, where Sg is plotted as a function of AT = — T for two fixed typical lattice temperatures T = 4 and 40 K 
of the two regimes. It is seen that in LLTR, the hot-electron effect first increases Sg very rapidly as rises from T 
to about several hundred degrees Kelvin (104 K in the inset for T = 4 K) and then gradually reduces it if goes on 
increasing. Nevertheless, even though is as high as 1004 K as shown in the inset of Fig. 6, Sg is still so large at 
low lattice temperature that it can not be neglected as suggested by Xing et al£3. However in HLTR, Sg decreases 
slowly as Tg increases. The results for Sg in regard to hot electron effects are somewhat different at higher density for 
the sample of Fig. 7 with = 3.6 x 10^^ m~^. In this case, Sg is plotted versus T as a solid curve for = T + 100 
K and as a dashed one for = T + 10 K, with all the other material parameters being the same as those used in 
Fig. 4. For comparison we also plot Sg in the linear transport condition as a dotted curve in the same figure, which 
has been obtained in Fig. 4. Again Sg is markedly enhanced in LLTR for T < 10 K. However, for Te = T + 100 K 
(solid curve), Sg is now very close to the plot for Tg = T -f 10 K (dashed curve) in LLTR and Sg seems to be almost 
unchanged in HLTR. This is further illustrated in the inset of Fig. 7, where Sg is plotted as a function of AT — T^ — T 
for three fixed lattice temperatures T = 3, 6, and 15 K. It is seen that Sg increases to "saturation" steeply in LLTR 
as Te increases and then rises very slowly as electron temperature goes up to several hundred degrees Kelvin. As Tg 
further increases, Sg begins to decrease slowly. Notwithstanding the fact that Tg is as high as T -I- 1000 K and Sg is 
slowly decreasing, the absolute value of Sg is still much higher than that in the linear-transport limit (AT — 0). In 
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HLTR, Sg is almost unchanged as Tg rises to several hundred degrees Kelvin before it begins to decrease. These new 
features result from the high density of the electron gas. 

It is well known that the phonon-drag component of thermoelectric power is very important in the linear transport 
limit and is usually so large that the difEiision contribution would not be discernible unless measurements were extended 
to very high or very low temperaturcsll3. Both the theoretical calculations of Lcitil and Xing et aLll3 show that the 
hot-electron effect does not change the absolute value of the diffusion component of thermoelectric power significantly. 
However, our results show that it is necessary to include the phonon-drag component Sg under hot-electron transport 
conditions, provided that the lattice temperature is not too high, and the electron drift velocity is smaller than the 
sound velocities, to accurately calculate the thermoelectric power. 

In summary, we have determined the role of phonon-drag in thermoelectric power under hot electron transport 
conditions using balance equation theory for weakly nonuniform systems. Having introduced wavevector- and phonon- 
mode-dependent relaxation times due to phonon-phonon interactions, the formula obtained can be used not only at 
low lattice temperature (where the phonon relaxation time is determined by boundary scattering), but also at high 
temperatures where phonon-phonon scatterings are dominant in determining the phonon relaxation time. Our results 
compare well with experimental data in the linear transport .-limit for a wide range of electron densities and are in 
good correspondence with the earlier theoretical predictionQM of the Boltzmann equation, appropriately modified. 
The role of hot-electron effects on the phonon-drag contribution to thermoelectric power Sg has been clarified, and the 
importance of the contribution of Sg to thermoelectric power under hot-electron transport conditions was stressed. 
In fact, the sign of the total thermoelectric power S in the presence of a sufficiently strong external electric field 
depends on competition between the negative Sd and positive Sg values for high impurity-density materials. For such 
materials, at least for the low lattice temperatures, the negative S^ value may be entirely counteracted by Sg, so that 
the total thermoelectric power remains positive. 
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FIGURES 



FIG. 1. Variation of thermoelectric power S* as a function of temperature T. The black dots are the experimental 
data of Sample 1 in Ref. ^. The soHd curve is the theoretical result obtained from Eqs. (31), (35) and (34) with Ax 
in Eq. (|l^) chosen to fit the experiment. The components due to transverse piezoelectric coupling Sgpt (short-dashed 
curve), longitudinal deformation potential coupling Sgi (dotted curve) and longitudinal piezoelectric interaction Sgpi 
(long-dashed curve), together with the diffusion component Sd (chain line), are plotted in the same figure. The 
parameters we employ are listed in the text. 

FIG. 2. Thermoelectric power S vs. temperature T. The black dots are the experimental data of Sample 2 in Ref. 
^. The solid curve is our theoretical prediction using the same values of A\ in the calculation as those obtained from 
Fig. 1. The components of S are represented by the other curves in the same way as in Fig. 1. 

FIG. 3. Thermoelectric power S is plotted as a function of temperature T. The black dots refer to the experimental 
data of Sample 3 in Ref. ^. The other curves represent components of S as in Fig. 1. 

FIG. 4. Comparison of thermoelectric power as given by the-itheoretical prediction with experimental data. The 
open circles are the recent experimental data of Fletcher et a/.ll3 . The curves represent S and its components as in 
Fig. 1. 

FIG. 5. Comparison of the Boltzmann equation predictions of Fletcher et (with the approximation explained 
in the text) with our balance equation predictions for S. In the figure the chain curve is the theoretical prediction of 
our balance equation formulation; the solid curve is obtained from our formulas (pT|), ( |35| ) and (^), however, with the 
approximation of Cantrell and ButcherQ installed as explained in the text, following Fletcher et al.; and the dotted 
curve is taken from Ref. |l^ as the theoretical prediction of Fletcher et al. 

FIG. 6. Phonon-drag component Sg vs. lattice temperature T for the sample N — 1.78 x 10^^ m^^. Solid curve: 
in the hot-electron transport condition ~ T + 200 K; dashed curve: = T + 100 K; chain curve: Tg = T + 10 K; 
and dotted curve: in the linear transport limit. Inset: Sg is plotted as a function of AT = T^ — T for two fixed lattice 
temperatures T = 4 K and 40 K. The material parameters are the same as those used in Fig. 1, as described in the 
text. 

FIG. 7. Phonon-drag component Sg vs. lattice temperature T for the sample N = 3.6 x 10^^ m~^. Solid curve: in 
the hot-electron transport condition Te = T -t- 100 K; dashed curve: = T -I- 10 K; and dotted curve: in the linear 
transport limit. Inset: Sg is plotted as a function of AT = Te — T for three fixed lattice temperatures T = 3 K, 6 K, 
and 15 K. The material parameters are the same as those used in Fig. 4 as described in the text. 



8 



